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The Gerasimov-Drell-Hearn integral Igdh(Q 2 )> an d its relation to polarized nu- 
cleon structure functions, is discussed from the lattice perspective. Of particular 
interest is the variation of Igdh(Q 2 ) with Q 2 , and what it may teach us about 
the origin and magnitude of higher- twist contributions. 



1. Introduction 

The Gerasimov-Drell-Hearn (GDH) integral, which is written as 

dv r 



GDH 



(Q 2 



87r 2 a 



a-{v, Q z ) 
1 



■»'<) 



Q 2 jo 

167T 2 a 



dx ■ 



dx 



5i (x, Q 2 ) 



a^{v, Q 2 ) 



(1) 



7 2 S2(*,Q 2 ) 



where v a — m x + (m 2 + Q 2 )/2mjv, xo = Q 2 /2mN^o and 7 2 = Am 2 N x 2 /Q 2 , 
connects the GDH sum rule at Q 2 — to the Bjorken and Ellis-Jaffe sum 
rules at large values of Q 2 . The spin asymmetry A\ is known over a large 
kinematical region for proton, deuterium and helium targets, which allows 
to compute /gdh(Q 2 ) down to Q 2 w 1 GeV 2 , separately for the proton and 
the neutron. 
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The GDH integral is of phenomenological interest for several reasons. 
It involves both polarized structure functions of the nuclcon, g\ and g 2l 
and thus tests the spin structure of the proton and the neutron. Further- 
more, the GDH integral provides a link between the nucleon state at high 
and at low resolution, allowing us to study the transition from an assembly 
of quasi- free partons to strongly coupled quarks and gluons. In particu- 
lar, we hope to learn about the structure and magnitude of higher-twist 
contributions. This requires, however, that higher-twist contributions set 
in gradually and before Q 2 reaches « 1 GeV 2 , that is before the opera- 
tor product expansion breaks down. In order to match the predictions for 
^gdh(O), the GDH sum rule, with the Bjorken and Ellis- Jaffe sum rules, a 
strong variation of /gdh(Q 2 ) with increasing Q 2 is anticipated. 

Lattice QCD is in the position to address these questions. In this talk 
I shall confront measurements of Igdh(<3 2 ) with recent lattice results. 



2. Polarized Structure Functions 

Let me recapitulate what we know about the nucleon's polarized structure 
functions g\ and g2, which enter in (1), first. 

A direct theoretical calculation of structure functions is not possible. 
Using the operator product expansion, we may relate moments of structure 
functions in a twist or Taylor expansion in 1/Q 2 , 

f 1 1 

2 jf dxx n gi (x,Q 2 ) = -e hn (Q 2 /^ 2 ,g(^ 2 ))a n (^+0(l/Q 2 ), (2) 



2 / dxx n g 2 (x,Q 2 ) = n [e 2 , n (Q 2 /[i 2 ,g(v 2 ))d n (fi) 
Jo + -v 

- e hn (Q 2 /n 2 ,g(n 2 ))a n (v)] +0(1/Q 2 ), 



(3) 



to certain matrix elements of local operators 

(P, s\°{a^-^ n }\^ «) = ~p[ a ^ WPv* ■ ■ -P»n + traces] , 

+ ■ ■ ■ — traces] , 

where 



^Mi-Mn = ( 2/ 97<t75jD ^ '" jD ^<7 -traces. (5) 
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In parton model language 

al = 2 f 1 dxx n [qi(x,^)- qi (x,fi 2 )}=2A n q, (6) 
Jo 

V „ ' 

A<z(x, M 2 ) 

in particular Oq = 2Au, = 2Ad, while c% has twist three and no parton 
model interpretation. 

In the following I will restrict myself to nonsinglet and valence quark 
distributions due to lack of space. These quantities show little difference 
between quenched and full QCD calculations, so that I can further restrict 
myself to quenched results. 

Let us first look at the structure function g±. In Table 1 I compare the 
lattice results for the lower moments of Aq(x, Q 2 ) 1 with the corresponding 
phenomenological (experimental) numbers 2 . The quoted result for Au — 
Ac? = gA has been taken from a recent, 'proper' extrapolation to the chiral 
limit 3 , shown in Fig. 1. By and large we find good agreement. 

Let us next look at the structure function g 2 - This differs from gi 
by twist-three contributions. From (3) we readily see that g 2 fulfills the 
Burkhardt-Cottingham sum rule 

/ dxg 2 (x,Q 2 ) = 0. (7) 
Jo 



Moment 


Lattice 1 


Experiment 2 


Au v 


0.889(29) 


0.926(71) 


Ad v 


-0.236(27) 


-0.341(123) 


A x u v 


0.198(8) 


0.163(14) 


AH V 


-0.048(3) 


-0.047(21) 


A 2 u v 


0.041(9) 


0.055(6) 


A 2 d v 


-0.028(3) 


-0.015(9) 


Au - Ad 


1.25(7) 


1.267(142) 


A 1 !! - A 1 d 


0.246(9) 


0.210(25) 


A 2 u - A 2 d 


0.069(9) 


0.070(11) 



Table 1. Comparison of lattice and experimental values of the lower moments of 
A.q(x, Q 2 ), defined in (6), in the MS scheme at Q 2 = 4 GeV 2 . The subscript v refers to 
valence quarks. 
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Figure 1. The axial coupling as a function of together with a fit from chiral 
perturbation theory. 



The first nontrivial moment of the twist-three contribution, d\, can be 
related to the tensor charge Sq of the nucleon 4 , 

d{ {s^p v - s„p M ) = (p, s|g(7 M 7 5 D v - D^^ 5 jq\p, s) 
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where m q is the mass of the quark. Thus we have 

d?(Q 2 ) = ^«y ? (Q a ), 

m N 



(9) 



which vanishes in the chiral limit (m q — > 0). The second moment, d\, has 
been computed on the lattice 5 . The result is shown in Fig. 2, separately 
for the proton and the neutron, and found to be in good agreement with 
experiment 6 . From (2) and (3) we obtain 



f 1 2 f 1 1 

/ dxx 2 g 2 (x, Q 2 ) + - / dxx 2 g 1 (x,Q 2 ) = -d 2 
Jo ^ Jo o 



(10) 



Given the fact that c? 2 is small, and d\ even vanishes in the chiral limit, we 
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Figure 2. The twist-three contribution <fe(M) at £t 2 = 5 GeV 2 for proton (O) and neu- 
tron (#) plotted against the lattice spacing a, together with the experimental numbers 
for proton (□) and neutron (■). 



derive that the Wandzura-Wilczek relation 7 

/ dxx n g 2 (x,Q 2 ) = I dxx n 9l (x,Q 2 ), (11) 

Jo n+lj 

g 2 (x,Q 2 )= [ —gi(y,Q 2 )-gi(x,Q 2 ) (12) 
j x y 

holds to better than 0(5%), except perhaps at very large x, which we are 
not interested in here. It is needless to say that g 2 in (12) satisfies the 
Burkhardt-Cottingham sum rule as well. The structure function gi(x, Q 2 ) 
is obtained from the parton distributions Aq(x, Q 2 ) by 

9i^Q 2 ) = \Y, e l f-e^Q^Aql-M 2 ) (13) 

with 

/ dyy n e 1 (y 7 Q 2 )^e hn (l 7 g(Q 2 )). (14) 
Jo 



3. Higher-Twist Contributions 

Contributions of twist four (and higher) have been studied on the lattice, 
either through calculations of appropriate nucleon matrix elements 8 , or 
by evaluating the operator product expansion directly on the lattice 9 . 
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Figure 3. Diagrams involving four-fcrmion operators. 

Higher-twist contributions are generally found to be small. Four-fermion 
operators, for example, schematically drawn in Fig. 3, which account for 
diquark effects in the nucleon, contribute 

/ dx F 2 (4) Q 2 )\ 27 J=1 = -0.0005(5) a s (Q 2 ) ^ . (15) 
Jo *°t 

The reason for quoting the flavor 27-plet structure function here, rather 

than the nucleon (octet) one, is that the underlying four-fermion operator 

does not mix with operators of lower dimension, which makes it a clean 

prediction. 



4. The GDH Integral 

We are ready now to address the GDH integral. Following our previous dis- 
cussion, we may assume that higher-twist contributions are small, and that 
the nucleon's second polarized structure function gi is well approximated 
by the Wandzura-Wilczek form (12). Moreover, given the fact that the 
lattice predictions for a%, compiled in Table 1, are in good agreement with 
the phenomenological numbers quoted, we may base our further discussion 
on the parameterization of gi(x, Q 2 ) given in Ref. 2. 

As already mentioned, I will restrict myself to the nonsinglet GDH in- 
tegral, which corresponds to proton minus neutron target. In Fig. 4 I show 
recent results from the Hermes collaboration 10 . I compare this result with 
the theoretical predictions. The solid line represents the full integral, as 
given by the bottom line of (1), including gi and g 2 , while the dashed line 
represents the integral 

167r 2 a f 1 , , ~ 167r 2 a , N 

^ dxg^Q 2 )^-^-^. (16) 
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Figure 4. The GDH integral: Experiment versus theory. The inner error bars are 
statistical, the outer ones statistical plus systematic. 
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Figure 5. The GDH integral: The connection to the GDH sum rule. The left curve is 
the prediction of chiral perturbation theory. 
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Equation (16) corresponds to the limiting case xq = 1 and 7 = 0. The 
dashed curve can hardly be distinguished from the solid curve over the 
whole kinematical range, 1 < Q 2 < 00, which tells us that the (nonsinglet) 
GDH integral is insensitive to g 2 and merely tests the Bjorken sum rule 



Jo 



dxgt(x,Q 2 ) = ^g A 
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(17) 



At larger values of Q 2 the data fall below the curve. This may be due 
to the fact that the x range covered shrinks to 0.2 < x < 0.8 in the highest- 
Q 2 bin. In order to draw quantitative conclusions, one certainly would 
need more precise data. So far we can only say that the experimental data 
are consistent with leading-twist parton distributions all the way down to 
Q 2 w 1 GeV 2 , and with the Bjorken sum rule in particular. 

In Fig. 5 I show the same figure on an enhanced scale, together with 
the predictions of chiral perturbation theory 11 and the predicted value 
of the GDH sume rule. Both curves appear to meet in a sharp peak at 
Q 2 « 0.4 GeV 2 . This value lies much below the range of validity of the 
operator product expansion. Likewise it lies beyond the applicability of 
chiral perturbation theory, so that no firm statement about the transition 
from large Q 2 to the resonance region can be made. 



5. Conclusions 

To learn anything new, if possible at all, from the GDH integral, we need 
better experimental data. At present the GDH integral does not teach 
us anything quantitative about higher-twist contributions. Lattice calcula- 
tions, on the other hand, indicate that the GDH integral is well represented 
by the Bjorken sum rule down to Q 2 « 1 GeV 2 . More precise lattice data 
on moments of gi and 52, including sea quark effects, will become available 
soon 12 . In order to match scaling and resonance region, it appears that 
one will have to resort to model building. 



Acknowledgement 

I like to thank Johannes Blumlein, Helmut Bottcher, Meinulf Gockeler, 
Roger Horsley, Dirk Pleiter and Paul Rakow for discussions and collabora- 
tion. 



9 



References 

1. M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz 
and A. Schiller, Phys. Rev. D53, 2317 (1996); M. Gockeler, R. Horsley, L. 
Mankiewicz, H. Perlt, P. Rakow, G. Schierholz and A. Schiller, Phys. Lett. 
B414, 340 (1997); S. Capitani, M. Gockeler, R. Horsley, H. Perlt, D. Petters, 
D. Pleiter, P.E.L. Rakow, G. Schierholz, A. Schiller and P. Stephenson, Nucl. 
Phys. Proc. Suppl. 79, 548 (1999); M. Gockeler, R. Horsley, W. Kurzinger, 
H. Oelrich, D. Pleiter, P.E.L. Rakow, A. Schafer and G. Schierholz, Phys. 
Rev. D63, 074506 (2001); M. Gockeler et al, in preparation. 

2. J. Bliimlein and H. Bottcher, Nucl. Phys. B636, 225 (2002). 

3. T.R. Hemmert, M. Procura and W. Weise, Phys. Rev. D68, 075009 (2003). 

4. This derivation is due to Paul Rakow. Equation (9) is still to be verified 
numerically. 

5. Last reference in Ref. 1. 

6. PL. Anthony et al. [E155 Collaboration], Phys. Lett. B458, 529 (1999). 

7. S. Wandzura and F. Wilczek, Phys. Lett. B72, 195 (1977). 

8. M. Gockeler, R. Horsley, B. Klaus, D. Pleiter, P.E.L. Rakow, S. Schaefer, A. 
Schafer and G. Schierholz, Nucl. Phys. B623, 287 (2002). 

9. S. Capitani, M. Gockeler, R. Horsley, H. Oelrich, D. Petters, P. Rakow and 
G. Schierholz, Nucl. Phys. Proc. Suppl. 73, 288 (1999). 

10. A. Airapetian et al. [Hermes Collaboration], Eur. Phys. J. C26, 527 (2003). 

11. V. Bernard, T.R. Hemmert and U.G. Meissner, Phys. Lett. B545, 105 (2002). 

12. M. Gockeler et al., in preparation. 



